THE Dfc-HEAT EQUATION ON QUADRIC MANIFOLDS 

ALBERT BOGGESS AND ANDREW RAICH 

Abstract. In this article, we give an explicit calculation of the partial Fourier transform 
of the Dfe-heat equation on quadric submanifolds of M C C™ x C m . As a consequence, we 
Q\ ' can also compute the heat kernel associated to the weighted 9-equation in C n when the 

weight is given by exp(— <f>(z, z) ■ A) where (f> : C™ x C™ — * C m is a quadratic, sesquilinear 
form and A € R m . Our method involves the representation theory of the Lie group M and 
the group Fourier transform. 
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1. Introduction 



> 

The purpose of this article is to present an explicit calculation of the Fourier transform of 
the fundamental solution of the LJ&-heat equation on quadric submanifolds M C C n x C m . 
A quadric submanifold can be thought of as a generalization of the Heisenberg group - 
it is a Lie group with a known representation theory |PR03j . and the technique of using 
Hermite functions to compute the heat kernel, as done in |Hul76j IBR 09J and elsewhere, can 
be extended to work in this situation as well. 

A consequence of our fundamental solution computation is that we can explicitly com- 
pute the heat kernel associated to weighted <9-problem in C n when the weight is given by 
exp(— <p(z, z) ■ A) where <fi : C n x C n — > C m is a quadratic, sesquilinear form and A G R m . 
This computation partially generalizes the results in [BR09J. When m = 1 and the weight is 
given by exp(rP(zx, . . . , z n )) where r G M, P(zi, . . . ,z n ) = ^2^ = xPj(zj), and pj are subhar- 
monic, nonharmonic polynomials, Raich |Rai06al IRai06bt IRai07l IRaij has estimated the heat 
kernel associated to the weighted (9-problem. If, in addition, n = 1, the weighted 9-problem 
and explicit construction of Bergman and Szego kernels have been studied by a number of 
authors in different contexts, e.g., |Chr9H IHasM IHas95l IHas98l lFS9~Tl IBer92j . We also note 
that quadric manifolds are related to if-type groups on which Yang and Zhu have computed 
the heat kernel for the sub-Laplacian [YZ08j . 

The remainder of the paper is organized as follows: in Section [21 we define our terms and 
state our main results. Section [3] provides the necessary background from representation 
theory. In Sections H] and 0, we apply the representation theory to the heat kernels and 
prove the main results. 
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2. QUADRIC SUBMANIFOLDS AND THE Dfe-HEAT EQUATION. 

2.1. Quadric submanifolds. Let M be the the quadric submanifold in C n x C m defined 
by 

M= {(z,w) e C n x C m ; Imw = 4>(z,z)} 

where <fi : C n x C™ i— ► C m is a sesquilinear form (i.e. 0(z, z') = <p(z', z)). For emphasis, we 
sometimes write to denote the dependence of M on the quadratic function 0. Note that 
M_ ( p is biholomorphic to by the change of variables (z,w) i— > (z, — w). 
For A G M m , let 

A (z,/) = 0(z,/) • A 

where • is the ordinary dot product (without conjugation). Observe that <p x (z,z') is a 
sesquilinear scalar-valued form with an associated Hermitian matrix. Let v x , . . . , v x be an 
orthonormal basis for C n with 

where //j(A) = /i A are the eigenvalues of the matrix associated with <j) x . 

2.2. Lie Group Structure. By projecting M C C n x C m onto G = C n x M m , the Lie group 
structure of M is isomorphic to the following group structure on G: 

gg'= (z, t) (z', t') = (z + z', t + t' + 2 Im </>(z, /)) . 

Note that (0, 0) is the identity in this group structure and that the inverse of (z, t) is (—z, —t). 

The right invariant vector fields are given as follows: let g £ G; if X is a vector field, 
then we denote its value at g by X(g) as a element of the tangent space of M at g. Define 
R g : G i— > G by R g (g') = g'g; then the right invariant vector fields, X(g), are obtained by 
pushing forward the vectors in the tangent space at the origin via the differential of the map 
R g . In particular, a vector field X is right invariant if and only if X(g) = (R g )*{X(0)}, 
where (R g )* denotes the push forward operator. Let v be a vector in C n ~ M. 2n which can 
be identified with the tangent space of M at the origin. Let d v be the real vector field given 
by the directional derivative in the direction of v. Then the right invariant vector field at an 
arbitrary g — (z,w) £ M corresponding to v is given by 

X v (g) = d v + 2 Im <f>(v, z) ■ D t = d v — 2 lmcj)(z, v) ■ D t 

where D t = (d tl , . . . , d tm ), (see Section 1 in Peloso/Ricci |PR03j ). Let Jv be the vector in 
R 2n which corresponds to iv in C n (where i = a/— T). The CR structure on G is then spanned 
by vectors of the form: 

Z v {g) = (1/2) (X v - iX Jv ) = (1/2) (d v - id Jv ) - icj){z~v) ■ D t 

and 

Z v (g) = (1/2)(X V + iX Jv ) = {l/2){d v + id Jv ) + i<f>(z, v) ■ D t . 

Also, 

[X V ,X V ,\ =Um<f>(v', v) ■ D t , [Z v , Z v ,] = 

and 

[Z v , Z v >] = 0, [Z v , Z v ] = 2i(f)(v, v') ■ D t . 
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We often drop the g in the vector field notation. The vector field definition of the Levi form 
of M is the map v i— > proj([Z„, Z v ]), where proj stands for the projection onto the totally 
real part of the tangent space of M at the origin (i.e. the t-axis). From the above equation, 
clearly the Levi form of M can be identified with the map v i— > <p(v, v), as mentioned at the 
beginning of this section. 



Recall that for any A G 



diagonalizes 



the set of vectors v*, 
A. For A G 



, v„ is an orthonormal basis which 



l , define the function z/(A) by 
u(X) = rank(0 A ). 
The function z/(A) satisfies < u(X) < n and as in |PR03| . 



{A G M m : u(X) = maxi/(A)} 

AeM m 



that carries full measure, i.e., | 
■ ^ We also write z = YTj=\ [ 



\Q\ =0. We identify x with 

A + iyfyvj for z = x + iy G C n . 



is a Zariski-open set f2 C R 

. . . , x A ) and ?/ with (?/ A , . . . , y^) We also write z = }_^'j =1 (x 
Additionally, we let z' = (z±, . . . , z\ X s), z" = (^( A ) +1 , • • • , z£) and similarly for x and y. 

Since the right invariant vector fields corresponding to <fi are equal to the left invariant 
vector fields corresponding to — and M_0 is biholomorphic to M^, any analysis involving 
right invariant vector fields yields corresponding information about the left invariant vector 
fields and vice versa. 



2.3. □& Calculations. Let v±,...,v n be any orthonormal basis for C n . Let Xj = X Vj , 
Yj = Xj Vj , and let Zj = (1/2) (Xj — iYj), Zj = (1/2) (Xj + iYj) be the right invariant vector 
fields defined above (which are also the left invariant vector fields for the group structure with 
4> replaced by — 0). Also let dzj and dzj be the dual basis. A (0, g)-form can be expressed 
as Y^Kelg ^ K wnere 1 q — {K = (ki, . . . , k q ) : 1 < ki < ■ ■ ■ < k q < n}. Proposition 2.1 
m [PR03] states that 

B b (J2 ( t ) Kdz K )= D LK^ K dz L 



Kei„ 



K,Lei„ 



where 



with the sub-Laplacian on G 



n LK = -5 LK C + M LK 



£ = (1/2) J2 z kZ k + Z k Z h 



and 



M, 



LK 



k=l 

^^[Zk, Zk] — ^^[Zk, Zk\ 

\k£K k^K 

e(K,L)[Z k ,Zi\ 




if K = L 

it\KnL\ = q-l 
otherwise 



Here, e(K,L) is (— l) d where d is the number of elements in K H L between the unique 
element k G K — L and the unique element / G L — K. The above theorem is stated and 
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proved in [PR03J for the left-invariant vector fields. If right invariant vector fields are used, 
then the above theorem provides a formula for associated to M_ t j > . 

For later, we record the diagonal part of i.e., D LL . Using ([T]) with L = K and the 
above formulas for Z k , we obtain 

i ( n ) n 

D LL = -- A + 2 Im <^> v *) d *u r D t ~ W^> Vk ^ ' Dt ) M*' Vk ^ ' Dt ) 

lfc=l J k=l 

+» y] 0(«fc, «*) • a - «*) ■ a ) (2) 

where A is the usual Laplacian in the ^-coordinates. For example, in the classic case of the 
Heisenberg group, (f>(z, z) = \z\ 2 , and Z k = d Zk — i~zj:d t , and D b is a diagonal operator (since 
[Z k , Zj\ = when k ^ I). The above formula for then gives the coefficient of acting 
on forms of the type <p L (z)dz L . 

2.4. The Dfj-heat equation and the Fourier transform. The heat equation defined on 
(0, q) -forms on M is the initial value problem on s G (0, oo) and (z, t) G M given by 

^ + U b p = in (0, oo) x M 

OS 

p(s = 0, z, t) = 6 Q (z, t) on {s = 0} x M 

Here, s is the time variable and t G M m is a spatial variable. Although we cannot find a closed 
form for p(s, z, t), we can find the partial Fourier transform of p(s, z, t) in the t- variables. 
Given a variable t G R, the (partial) Fourier transform in t is given by 

V •/ K 

If / is a function of several variables f(t\, . . . ,t k ) and, for example, we take the partial 
Fourier transform in ti, we use the notation f(r, t 2 , • • • , t k ). 

As we will see below, to compute the partial Fourier transform of p(s, z, t) = p s (z, t), it is 
enough to solve the (Fourier transform of the) DLL-heat equation 

|^ + D LL p = in (0, oo) x M ^ 
p s =o{z, t) = 5 {z, t) on {s = 0} x M 

We start by computing the partial Fourier transform in t of denoted 0^ L . We start 
with a reexamination of d2J). By taking the partial Fourier transform in t of the formula for 
□ll, the effect is to replace D t with zA. If we write z = Ylk=i z k v k, then 

n n n 

Im | ^2 ^(^ v k)d Zk | • zA = (f)(v k , v k ) ■ iX Im {z k d Zk } = ^ ip k Im {z k d Zk } 
fe=i fe=i fe=i 
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and 



^2 (0(z> v k) ■ iX) ((p(z, v k ) ■ iX) = ^ {zk<t>{vk, v k ) ■ iX) (z k (p(v k , v k ) ■ iX) = - ^{^ x k ) 2 \z k \ 



k=l k=l k=l 

Consequently, (J2J) transforms to 



□L = -|A + 2i^^Im {z k d Zk } + $>£) 2 N 2 - - ) W 

fc=l k=l \fceL fc0L / 

We employ the following notation: for 1 < j < u(X), define (L) = = sgn(/x* ), if j G L 
and 6j = —sgn(fij) if j $ L. 

Our main result is the following. 

Theorem 1. For any X G lR m , the partial Fourier transform to the fundamental solution to 
the D LL -heat equation satisfies the heat equation 

^ + n x LL p = in (0, oo) x C n 

p(a = 0, z, A) = (27r)- m / 2 5 (2) on {s = 0} x C n 

and given by 

p(*,ar,2/,A) = [ -=%, e-^^ TT w ^ , 

^ ' y 7 S n-^A) 11 smh s/i A 

Note that /z* and coth(s/i A ) are real-valued and are odd in /i A , so putting absolute values 
around the fij would not change the result. Therefore, there is Gaussian decay in (x 2 + yj) 
for all j when A G M m . Theorem [T] generalizes Theorem 1.2 in [BR09] in the case that r G K. 
and 7 = n — 2q. 

We now cast the heat equation in terms of a weighted <9-problem in C n . Recall that 
Zj = ^7 + i(j){z,v) ■ D t . If we denote a superscript A for the partial Fourier transform in i, 
then 

Z 3 ^z] = -^- (P(z, Vj ) ■ X = e *(*.*)-*A e -*(*.*)-A. 

3 3 

From the computation of Zj, the tangential Cauchy-Riemann operator is defined on 
(0, g)-forms on G by 



w*)= E , E<$!z j fj(z)dz 

where 



K ^ q +i j=i 

JGIn 



e jJ 



-l) a if {j} U J = K and a is the sign of the permutation taking {j} U J to K 
otherwise 
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This means that if g is a (0, g)-form in C™ and we treat A as a parameter, then the partial 
Fourier transform in t of db, denoted by d b is given by 

d X b g{z) = e#*'*)- A 0{e-*<*'*>' A 0} 

where d is the usual Cauchy-Riemann operator on C n . Since □& = <9b<9£ + d 6 d& where 9 6 is 
the L 2 -adjoint of d&, it follows that D b = d b (d b ) + (d b )*d b . Thus, solving for the D A -heat 
kernel also yields the heat kernel associated to the weighted <9-problem on C n with the weight 

e -<j>(z,z)-X 

Corollary 1. For any A G M m , the function 

H\s,z,z) = (27r) m/2 p s (^-5,A)e- 2iA - Im ^' i) 
satisfies the following: if 

H x {f}(s,z) = [ H x (s,z,z)f{z)dz, 

then H x {f} solves the initial value problem for the weighted heat equation: 

(d s + D x ){H x f} = zn(0,oo)xC" 
H x {f}(s = 0,z) = f(z) on {s = 0} x C n 

In particular, the component of H x (s, z, z) on dz L for L 6l, is 

nn~v{X)(c)\-n , „ .„ |2 u ^ O p S€ j \rf I ..A 

H£(s, z, z) = 2 (2 ,V e-^^ TT / ^ e-^^-)fe-%-l 2 e- 2 ^ Im ^). 

Note the formula for the heat kernel yields a standard Gaussian solution for the Euclidean 
heat kernel in the zero eigenvalue directions. Also, the disappearance of the (2ir)~ m ^ 2 owes 
to the fact that 5 (z,X) = (2^)- m ^5 (z). 

3. Representation Theory. 

3.1. Irreducible unitary representations. For z = x + iy E C n , t, A G M. m , and r\ E 
C n-KA) ? define ^ t) ; L 2 (R U ^) h-> L 2 (R U ^) by 

7r Ai7? (x,y,£)(/i)(e) = e i ( A -*+ 2R <- , '^)) e - 2 ^;a ) ^ 2 / A fe+^) /l(e + 2x /) 

for h G L 2 (lH A )) (so £ G M^ A )). Note that if r] = ( + iq, then Re(z" • fj)) = x" ■ ( + y" ■ q. 

The map ir\ >rj (x,y,t) is unitary on L 2 (R. U ^). Also, n is a representation for G, which 
means that for each A G f2, tt^,^ is a group homomorphism from G to the group of unitary 
operators on L 2 (R"( A )). Verifying that it\ jri is a representation is done in [PR03J. 

If A is a right-invariant vector field, then A gets "transformed" via 7T\ jV to an operator 
on L 2 (W^) denoted by T = d-n x ^(X). This means that 

x {n\M} = To7r x , v (g) (5) 

as operators on L 2 (M. U ^). It is usually easy to identify T by seeing what happens at g = 
and using the right invariance of A to show that the above equation holds for all g E G. 
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To clarify, let R g {g') = g'g and recall that the vector field X at the point g is given by 
X{g) = {R g ),{X{0)}. If X{vr Afl }(0) = T o 7r AjJ? (0), then we have 

X{n Xiri }(g) = (R 9 UX(0)}^xM} 
= X(g f = 0){n XtV (R g (g'))} 

= X(g' = 0){7r Ajr? ((/)7r A)r? (g)} since 7r is a homomorphism 

= { X (g' = °) 7r A,r)(5' / )} ° 
= T O 7T Afl (5() 

where the last equation uses the relationship of X(g) and it at g = 0. 
A similar computation shows that if X is left invariant, then 

X e {n Xtri }(g) = n Kn (g)oT 

as operators on L 2 (M^ A )). Note that the order of T and 7r A)J? is reversed from We will 
not dwell on this point as we prefer the use of right-invariant vector fields. The relationship 
X{nx tV }(g) = T o ir\ )V (g) is often expressed using the shorthand: dir\ tV (X) = T. 
From earlier, we have the right invariant vector fields 

Xj = d v x - 2 Im <f)(z, w A ) ■ D t 

and 

Yj = d Jv x - 2 Im (f>{z, iv}) ■ D t = d Jv x + 2 Re <f>{z, v$) ■ D t . 

where J is the usual complex structure map on M. 2n = C n . In view of (|3J), we have the 
following relations: for 

Y r k ^ f 2d ti°«xM l<j<u(\) 

Xji^vKg) = S o-/ n /, U1 . V ( 6 ) 
\2iQ o 7r\, v {g) v(\) + 1 <3 <n 

yr wx f-2i»jZj°*\M 1 < j < u(X) 

[2«sj 07^(0) v(\) + l<]<n 

dt k {^X,r t }(g) = ^k°TT\,ri(g) I < k < 1TL (8) 

as operators on L 2 (R^ A )). In the second equation, £j is thought of as a multiplication 
operator on L 2 (M^ A )), i.e. /(£) i— > Equations (jSJ) and ([7]) are easily shown to hold at 

the origin since Xj(0) = d x . and 5^(0) = d yj , and the right invariance forces these equations 
to hold at all g G G. 

Now we compute the "transform" of D LK in the coordinates (z A , . . . , z A ). Note that 

d-K X>v [z j ,z l } = { i > = ; 

This follows from (jSJ) and the fact that the coordinates (z A , . . . , z A ) were chosen to diagonalize 
the form <p(z, z) ■ A. In view of ([TJ and flBj) - (|Sj), we have 

M« = { " A < + "" 2 + ^KfZ - S2S«MI « * = * (9) 
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We will also need to transform the adjoint of Due which is defined as 

/ D LK {f(z,t)}g{z,t)dxdydt = f(z,t)Df^{g(z,t)} dx dy dt 

J(z,t)eG J(z,t)eG 

(note: this is the "integration by parts" adjoint, not the L 2 adjoint, since there is no conju- 
gation). We have 

== *,,nS = ( " A < + i"i 2 + + e£S ^ * = f do) 

? [ it K L 

(just a sign change for the last term on the right). The subscript £ on Q^' v ' indicates that 
this is a differential operator in the £ variable (instead of the group variable g = (x,y,t)). 
Below, we assume L = K (otherwise the operator is zero) and that L, A, and r] are fixed. 
We drop the superscript LL when its use is unambiguous. In view of ([T]), and (jSj) through 
(JHJ), we have 

Ofl{n x M} = Q^on x M (11) 
as operators on L 2 (R. U ^). We return to this key equation later. 

3.2. Group Fourier Transform. For (z, t) G G, we express (z, t) = (x, y, t) = (a/, y', x", y", t) 

(x', y', z", t). The variable z" may be thought of as in C n - u(x) or R 2 ( n ~»W). 

For / : G I— > C, we define the group Fourier transform of / as the operator Tj ,v : 
L 2 (R"W) ^ L 2 (R"W) where for h G L 2 (IT (A) ), 



Tj' v {h}(0 = [ f(z,t)ir XiV (z,t)(h){Odxdydt 



(z=x+iy,t)GG 

f 

(z=x+iy,t)£G 



/(^,t)e i ( A -* +2Re ( 2 ''^)e- 2i ^ ) ^^fe + ^)/ i (£ + 2x') dxdydt. 



As before, the coordinates for x, y G R n relative to the basis v A , . . . , v x . Note that 

T^ v {h}(0 = (2irY 2n+m ~ u ^/ 2 [ f(x', 2fi oTf+V), ^2 V , ^A>(£ + 2x') dx'. 

We have written /i A o (£ + x') for (/i A (£i + a; A ), . . . , ^l^(C,u(x) + ^(a)))- We can a ^ so ex P ress 
^'"{/i} as 

T^{h}(0 

= ( 2vr )(^+— (A))/2 /" ^„ ^{/(z, y, t)e~ 21 ^ ^}(x', 2/io£, -2r?, -X)h^+2x') dx'. 

(12) 

In the above notation, ^v',y,t indicates the Fourier transform in the (x",y,t) variables only, 
whereas JF indicates the Fourier transform in all variables (except s). 
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In view of (llljl . we have 

QpiTpih)®} = [ f(z,t)nfi{n x>v (z,t)h(Z)}dxdydt 

J (z=x+iy,t)eG 

= [ n LL {f{z,t)}TT X>r ,{z,t)h{Z)dxdydt (13) 

4. The Heat Equation. 

4.1. The heat equation on M. Our goal is to find a formula for the fundamental solution 
to the heat equation Q. We know abstractly that p exists: D LL is self-adjoint and nonneg- 
ative, so e~ sULL is a well-defined, bounded linear operator on L 2 {G) with norm at most 1. It 
has an integral kernel by the Riesz Representation Theorem, and the existence of p follows 
easily. See |Rai06a] for details. 

Let us apply the group Fourier transform to p and recall that p s (z,t) = p(s,z,t). Define 
the operator U*>*(s) : L 2 (R U W) h+ L 2 (W^) by 

U^(s){h}(0 = T^{h}(0 = [ p s (z,t)ny v (z,t)h(Odzdt. (14) 

J(z,t)eG 

In view of ( I13p and the fact that p s (z,t) solves the heat equation, we have 

Q X { v {U x ' r >(s){h}(0} = [ n LL {p s (z,t)}ir x>TI {z,t)h(Odzdt 

J(z,t)eG 



= -9 S / p s (z,t)ir x , v (z,t)h(£,)dzdt 
U(z,t)eG 

= -d s {u^( S ){hm} 

Also 

U^is = 0){/i}(£) = T£"{h}(£) = h{0- 
Therefore, we conclude that U x,n (s) satisfies the following boundary value problem: 

Q x ' n {U x (s)} = -d s {U\s)} and U x (s = 0) = Id (15) 

where Id is the identity operator on L 2 (M. U ^). This is a Hermite equation similar to, though 
more complicated than, the one we solved in the Heisenberg group case [BR09]. So, our 
approach is to proceed as follows: 1) explicitly solve this Hermite equation, and then 2) 
recover the fundamental solution to the heat equation. 

As to the second task, we let a G M^ A ) be an arbitrary vector, and then define h a (£) = 

( 27r )-n-m/2 e -i£-a Let 

u x «(s,a,0 = U x «(s){h a }(0- 

The above definition needs explanation since h a G" L 2 (IR ly ( A '). For each fixed s > 0, p s G 
L 2 (G) and we can approximate p s by p 5 s G L 1 R L 2 {G) (e.g., by multiplying p s with an ap- 
propriate test function). Then, as we see below, we can define U^ ,v (s){h a }(^) = T X g V {h a }(l;) 

Ps 
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since in view of equation ( fl2l) . 

u^(s){K}(0 

= , 9 w A)/2 / F*»,y,t U{x, y, t)e~ 2 ^ } (*', 2/i o i, -2rj, -\) e -^'> a dx> 
= T [p s s (x, y, t)e~ 2i ^ fa* } (2a, 2p x o £, — 2»|, -A)e^ a 
By the definition of the Fourier transform in L 2 , 
^ [p s s (x, y, t)e~ 21 ^ fa'i } (2a, 2// o £, -2ij, -\)e** 

— ^ {p s (x, t)e- 2 * ^ | (2a , 2/i A o £, -2tj, -A) e -* a 

in L 2 (R. U ^) as 5 — > 0. Thus, u x,n (s,a, £) is well-defined. In the above computation, we view 
?7 = (C)<r) G R 2 ( ?1 ~ 1/ ( A )). Also, the motivation for the choice of h = h a is that it offers the 
"missing" exponential needed to relate the full Fourier transform of p s with u x ' v . Now it is 
just a matter of unraveling the equation 

u^(s,a,0=^{p s (x,y,t)e- 2 ^^ x i} (2a, 2p x o £, -2r/, -A)e"^ (16) 

for p s using the inverse Fourier transform. 

Before we go on, let us remark that had we used left invariant vector fields rather than 
right invariant ones, then the transformed operator, Q x,v would appear on the right of the 
group transform. That is to say, we would be trying to solve the following analogue of ({15]) 

ds{T^ s } = -T x Q^ and T x =o = Id 

where Q x ' v is a Hermite type differential operator similar to Q X ' V . Note the transform 
operator T x is now intertwined with the differential operators (i.e. d s is on the left side 
and Q^ v is on the right). Since the inversion formula for the group transform operator is 
complicated (see [PR03] ) . it would appear that this method of using left invariant vector 
fields is more difficult to unravel a formula for p. 

4.2. Weighted heat equation. Our objective is to compute p s (x' ,y',rj, A) by solving the 
weighted heat equation obtained by taking the partial Fourier transform in the t and (x", y")- 
variables. We obtain the D^-heat equation 

d s p s (x',y',r),\) = -U x L ;lp s (x',y',rj,X) 
p s= ,(x',y',%\) = (2rry m / 2 -^-^5 (x',y') 



From (J3J), we have 

„(A) u{\) u(\) 

j=i j=i j=i 
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In the following computation, we find a formula of p s (x',y',r}, A). Observe that pj x = — p x 
and ej x = — e x . (Note that vj x = v x , so we can continue to suppress the A superscript on 
Xj and yj.) We unravel ffTBl) to obtain (with a, b e R^ A )) 

Ps (x', y', % A) = e"*^ (e~i ^ ^ 6 ^^^( s , a , 6)) (a;', y') (17) 

where ■u A,r '(s, a, b) = u~ x, ~^ T1 (s, a/2, b/{2p~ x )) and b/(2p~ x ) is the vector quantity whose 
jth component is bj/(2p~ x ). As we shall see, the inverse Fourier transform in the a, and 
b variables will be relatively simple (using Gaussian integrals). In the next section, we use 
Hermite functions to solve for u x,ri on the "transform" side. Then we return to the above 
formula to compute p s (x', y', rj, A). 

5. Computing the heat kernels 
In this section, we prove Theorem [TJ and Corollary [TJ 

5.1. Hermite Functions. Our starting point is equation (|15|) . which we restate as Q x,v {U x,ri (s 
-d s {U x ' v (s)} where 

u(\) u(X) 

Q x , v = -A, + H 2 + J>:fe) 2 + E^l- 

3=1 3=1 

We use Hermite functions to solve this equation. For a nonnegative integer £, define 

Each ij}£ has unit L 2 -norm on the real line and satisfies the equation 

-^(x) + x 2 M%) = + l)^(^), 
see [Tha93| . (1.1.9). For A E R m \ {0}, define 

Each ip x .{^j) has unit L 2 -norm on R and hence ip x has unit L 2 -norm on R^ A ). An easy 
calculation shows that 

+ 0*&) a ){^&)} = + (18) 
For s > 0, we claim that [/^(s) : L 2 (R^ A )) \-> L 2 (W (X ^) as defined in (0 is given by 

f(A) oo 

where P £ A is the L 2 projection of a smooth function of polynomial growth in the variable 
£j onto the space spanned by ip x .(£,j), and where (S)ji A i is the tensor product (so that the 
output of U x ' v (s) is a function of £i, . . . , ^(x))- For shorthand, we write 

E x (s) = e-K^+^+^ll^l'. 
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We then have U x,ri (s) = e~ s ^ 2 (S^i ^4=0 -^4 ( s )-^y Using the product rule, we compute 



d s {U x «(s)} = -\ V \ 2 U^(s) + e-^J^ds E (g) £ E x k (s)P x 

3=1 V,=0 / £=i 4=0 

f(A) oo f(A) oo 

= -\v\ 2 U^(s) + e~^ 2 E -[(^ + !) + ^V [(2 ^ +1)+ * A ' S ^ <g) E 4(^4 

3=1 4=0 fc=i 4=0 

^(A) oo f(A) oo 

= -M 2 ^) + e- 1 - 12 E E (<U - (^0) 2 - ef 1^1) o e -[(^ +1 ) + ^]|^| S pA (g) £ ^a (s)p A 

i=i 4=0 fc=i 4=0 

where the last equality uses ( fTHI) . Since the differential operator on the right is independent 
of £j, we can factor it to the left of to obtain 

d s {U x «(s)} = -Q x «{U x «(s)}. 

Since the Hermite functions, t/>£, form an orthonormal basis for L 2 (IR), U x ' v (s = 0) is just 
the identity operator. Thus U x ' v (s) solves ffT5]) . 

As above, we apply U x ' r '(s) to the function h a (£) = (2ix)~ n ~ rn l 2 e~ l ^' a to obtain the funda- 
mental solution p s . We therefore obtain 

u x «(s,a,0 = U x «(s){h a (0} 

v{\) oo 

= (27r)— m / 2 e-^l 2 flj;4( 8 )^{e^} 

3=1 4=o 

Since ft, Q belongs to L°°(IR I/ ( A )) and not in L 2 (IR^ A )), the above sum converges a priori in the 
sense of tempered distributions (as opposed to L 2 convergence). Earlier, we argued that we 
can obtain U x,ri (s){h a } via a standard approximation argument, however, we will see below 
that the convergence is much stronger and the result is a smooth function in s, a, £. Each 
projection term on the right is 

P^e-^) = ( f e^^(lA*jr /2 &)l/'jr /4 d&) 
= (2 7 r)^( a ,/|^|^)^.(|^| 1 /2 e . ) 

= (27r) 1 / a H)^(o J -/|/*J| 1 / 2 )^(|^r /2 &) 

where the last equality uses a standard fact about Hermite functions that they equal their 
Fourier transforms up to a factor of (— i) 4 Substituting this expression on the right into 
the definition of u x,n (s, a, £), we obtain 

v(X) oo 
3=1 4=0 
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This function satisfies 

d s u^(s,a,0 = -Q^{u X ' v (s,a,0} 
u ^(a = o,a,0 = h a (0 = (2n)- n - m l 2 e- ia ^. 

In view of ( ITT)) , for computing p s (V, y',rf, A), we need to compute 

u*>*(s, a , 6) = u~ x '-^(s, a/2, b/{2^ x )) 

where bj (2/i~ A ) is the vector quantity whose jth component is bj/(2/j,J ). From the previous 
equality, and using that fij x = —fij, ej x = — e A , we have 

v(\) 

u x «( a ,a,b) = ( 2 7r)-5(«+("-^))) e - s ^ JJe-d-^l' 

3=1 

oo 

xE(^)' J ^^V2|/i J A r /2 )^,(^> J A | 1/ V2^ A )e-^^ |s . 

Let 



Do = e r , a, = - 
J 2|, 

Then 



C. = g-2|M?l« a . = ^ Q, = MO) 

bj 6 ' Qj %*|V*' Pj 2/i A • UJj 



2 "(A) . OO 

«^( S ,a,&) = (27r)-^" + -+(—W)) e -^ n^ 1_£ /2 E(-^i)^(^)^(/3i)- 
Using Mehler's formula ( |Tha93j . Lemma 1.1.1) for Hermite functions, we obtain 

„ . ,., _ - . ........ s^f- TT qC 1 -^)/ 2 1 -nT&J&i+W i+sT 



u A '"(s,a,&) = (2 7 r)-W 2 +")2^)/ 2 e-^n g i e 



The series for -u A,r? converges in C°° on the unit disk in C, and therefore the series for u x,v 
converges in C°° for s > 0, justifying many previous computations (which held a priori in 
the category of distributions). 

5.2. Finishing the proof of Theorem)]]. In view of Equation (|17|) . to determine p s (x', y', rf, A) 
we must compute 

Using (fl9l) . and simplifying we obtain 

e-iE^ a i b i/^u x ' v (s,a,b) = (27r)" (m/2+n) e~ s ^ II a,:„- ■ ■> m,,,,!,, 

j=i ^cosh(2|/i A |s) 

where 

t&nh(2\ u x \s) sinh 2 flu A ls) 
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After an exercise in computing Gaussian integrals, we obtain 



(27r) - (m /2 + n) e - S BL- JJ 



e j |P j 1 & y 3 3' 33 



After simplifying, 

"4> _ ..Xa in B i _o„a /_woi..ai^ /„2 ; R2 _ sinh( S |/iJ|) 



The previous expression becomes 
^(e^*WWfr*(8,a,bj) (x',y') 



,2 2p £ >j |s I/v a 
(27r) - (m / 2+ n) e - s ^ JJ ^ I/'/ 



fi sinh(s|^ A |) 

In view of ifTTj) . the fundamental solution p s (V, y', fj, A) to the weighted heat equation is ob- 
tained by multiplying this expression by YYj^i e~ 2lfJ, 3 x 3 y 3 which cancels the similar expression 
on the right side. We therefore obtain 



1 ,2 v{X) 2e e >j\ s \u x \ 
y', % A) = (2n)-^e-^ ]J - ^ 

A ^ smh(s|/i^|) 



" (A) 9^ a Im a I s I„ai 

l e -^(^/Bj)(xf+^)_ 

^ smh(s|^ A |J 
Note that the rightmost exponent can be rewritten as 

+ y]) = - f^^ tf + V]) = coth( M ^)(x| + yj). 

Consequently, 



p»^y^) = ^=ka) e " s 11 



^ ^coth(/i^s)(a;|+2/|)_ 



J sinh(s|/x A |) 
This completes the proof of Theorem [1] for A G Q. 

5.3. The proof of Corollary [TJ, In this subsection and the next, we show that the following 
kernel: 

H x (s, z } z) = (2ir) m/2 Ps (z - z, \) e -^™^) (20) 
= 2— ( A >(2vr)-" c _^ ^ 2e-»| / i A | c _ M , coth(M , s)|2i _ ii ,^_ 2tA , Im , ( ^ ) 



L sinh(s|/i A |) 



THE Db-HEAT EQUATION ON QUADRIC MANIFOLDS 15 

is the heat kernel for the weighted Here, z = x + iy and z = x + iy. Note that H is 
conjugate symmetric, i.e. H x (s,z,z) = H x (s, z,z). We will show that the heat kernel has 
the following properties: if / G L 2 (C n ), then 

H X {f}(s, x, y) = / H x (s, x, y, x, y)f(x, y) dx dy 

is the solution to the following boundary value problem for the heat equation: 
(d s + D x ){H x f} = 0, H x {f}(s = 0,x,y) = f(x,y). 

5.4. Group Convolution and Twisted Convolution. To motivate the above formula, 
we consider the fundamental solution to the (full) unweighted heat equation: p s (x,y,t). For 
a function / G L 2 (C n x R m ), and g = (z,t) eC"xl m define 

H{f }(s, g ) = ( Ps *f ){g) = [ps(g\g]- l )h(~g)d~g (21) 

J 9 

where * is the group convolution and g^] -1 is the group multiplication of g by the inverse 
of g. If X is a right invariant vector field, then 

XH{f }(s,g) = l{Xp s ){g\g]- l )U{g)dg. 

Since is comprised of right invariant vector fields and p s satisfies the Ob-heat equation, 
we therefore have 

(d s + D b ){H(f )} = 0. 
In addition, the following initial condition holds: 

H{f }(s = 0,g)= [p^gW^fM® = fo(9) 
J 9 

since p s= o(z,t) is the Dirac delta function centered at (z, t) = 0. 

Note that H x {f}(s, x, y) = (27i) m ^ 2 H{f}(s, x, y, A), which is the partial Fourier transform 
in the t variable of H{f}(s, x, y, t). We will now show the Fourier transform in the t- variable 
transforms the group convolution to a "twisted convolution" , which we now define. Suppose 
F and G are in L 2 (C n ), and A G M m . Following Stein |Ste93] . page 552, we let 

(F * A G)(z) = [ F(z- i)G{z)e- 2iX - lm ^~ z) dz. 

The arguments in [Ste93j . page 552 with (z, z) replaced by 2Im0(z, z) show the following: 
if F ,Gq G L 2 (C n x R m ), then 

(F * G Q )(z,X) = (2n) m / 2 (F (;X) * A G (-, X))(z). 

Now suppose / G L 2 (C n ) is given and let f (z,t) = (2n) m / 2 f(z)5 (t), so that f (z, X) = f(z). 
With H x given as in fTSOj) . we can take the partial Fourier transform in t of f[2T]) and use the 
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above relationship to obtain 

H x (f)(s,z) = [ H x (s,z,~z)f(~z)d~z 

= [ (27i) m / 2 p s (x -x,y-y, X)f(x, y) e -™-*"<K*# dx dy 

= (2vr) m / 2 (p s (-,A)* A / (,A))( 2 ) 

= (Ps*fo)(z,X) 
= H(f )(s,z,X) 

Since H(f ) satisfies the Da-heat equation, H(f )(s, z, A) = H x (f)(s, z) satisfies the weighted 
heat equation, i.e. 

(d s + n x b ){H\f)} = o. 

The initial condition H x (f)(s = 0, z) = f(z) is also satisfied because 

H x (f)(s = 0,z) = H(f o )(s = 0,z,\) 
= fo(z,X) 
= f(z). 
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